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ABSTRACT 


Water Ros;'’ur cos Systems are complos. Models are used 
for design and operation cf such ojmples systems. The models 
nay be physical cr mathematical. In the present study, a 
mathematical model is used for analysis of such systems. A 
set of simulation and optimization techniques used for analyzing 
the models of mu?Lt ire servo ir water resources system is explained. 
A resume of the development of synthetic streamf low generating 
techniques is given. Making use of an existing program for 
streamflow generation, synthetic streanflcws a.re generated for 
txfo neighbouring river basins. For optimization studies, the 
concept of representing the water resources system in a capaci- 
ta.ted network tom is made use cf. The .out~of -kilter algorithm 
which makes use of this concept is explained. An existing 
p-rogram, devol':o'’3d by the Texas Water Development Board, which 
makes use of the cut-c f-kilter al.go-rithm is used to study the 
performa.nce cf the miulti -re servo ir sy seem in one of the basins. 
The results of the simulation and optimization are presented 


and discussed 



CHAPTER 1 


IITRODUCTIOH 


1 . 1 General 

Water is one of the several resources without which a 
nation cannot satisfy the fundamental wants of its people or 
achieve the important national goals it sets for itself. 

Without water, life itself cannot he sustained. Our nation 
is blessed with a bovintiful supply of water, although it is 
not always in the right place at the right time. Because of 
the general abundance of water, it was taken for granted that 
water had no cost and there were no limits to its availability. 
But as demands come close to and in some regions even exceed 
supplies of water, it becomes necessary to seek ways to 
increase the efficiency in the use of water. The scope of 
this study is to optimally manage a complex water resources 
system involving a system of multiple reservoirs. The study 
would also reveal surpluses or shortages that would occur in 
the years of high or low flows compared with years of mean flows 

1 .2 Literature Review 

The application of computer oriented optimization technique 
to solve design and operation problems in the field of water 
resources development has received serious attention and support 
in the past few years. These techniques may be divided into two 
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main catagories i simulation and mathematical programming. 

Simulation has been effectively used in several investiga- 
tions to find the best feasible design for a complex multi- 
purpose multi-reservoir system (Maass, 1962). Hufschmidt and 
Piering (1966) have used the simulation technique to analyse 
Lehigh and Delaware water resources systems. 

The Lehigh Water-Resource System consisted of a system of 
six reservoirs out of which only five could be built at any 
time. Each reservoir was a multipurpose reservoir and 
could be used to provide (i) regulated flews for water supply 
or water quality improvement, (ii) recreation, (iii) flood 
protection, and (iv) hydro electric energy. The system also 
included a diversion channel. The simulation model of the Lehigh 
system had 42 major design variables. 

In these studies, a number of equations are formulated to 
describe the physical behaviour of the water resources system 
under investigation. Design variables (capacities of reservoirs, 
power plants, and irrigation canals) and operating rules 
(monthly release from reservoirs as a function of reservoir 
content and month of the year) are introduced into these 
equations as parameters. The problem is run through the computer 
for many combinations of the parameters. Physical feasibility 
of each design and the value of 1iie objective function are 
determined by the computer. The design with maximum value of 
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the objective function is chc son ss the optimuni design. 

These studios hnve shorn thp.t simulation technique perform 
s a t i s f a c t o r i ly . 

The introduction of the concept of reservoir zoning by 
Beard (1967), Fredrich and Beard (1972) was a significant 
development in the simulation technique. In this concept, 
every reservoir is divided into a number of storage zones 
and during simulation, all reservoirs are maintained in the 
same zone a,s far as possible. A generalized simulation program 
for the operation of a reservoir system for conservatiofi pur- 
poses such as water supply, navigation, recreation, low flow 
augumenta.ticn and hydro-electric power has boon deve3.oped by the 
Hydrologic Engineering Centre (1974). Visens and Schaake ( 1 980 ) 
have studied the Rio Colorado basin using the simulation 
technique . 

Another mpjor development in simulating complex water 
resources systems was the application of out~of -kilter algorithm 
to such systems (Ford and Fulkerson, 1'962) . The technique has 
been used by the Texas Water Development Board in a number of 
the ir program s ( T ¥DB - 1 , 1970, TWDB -2 , 1972). 

Ma.thernatical programming techniques are analytical methods 
with some theoretical assurance that the optimal solution will 
be reached given enough computer time. Some of the important 
works done using these to-chniques are reported in Chapter 3. In 
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aaciiticn tC' those 5 Harboe et al. (1970) used the dynamic 
progrejmning tc develop optimal po-licy for reservoir operation# 
To'rabi and Meba-sheri (l973) and lauze et al. (l980) have also 
used dj'naaic prugraaming for optimizaticn studies. Windsor and 
V.T. uhou (1972) used integer progranEiing end seperable 
programmirig tc develop a multi -^reservoir optiaiza,tion model. 
Himmc-lblau (1974) and Bayer (1974) used non-linear programming 
in their studies. 


In the present study, a monthly simulation program 
HEC‘-4 developed by Leo R. Beard (1972) of the Hydrologic 
Enginoering Centre, is used to- generate synthetic streamflows. 
.4n optimizoti: n model AL-3 , which makes use of the out-of- 
kiltor alg- rithj'n devoloped by the Texas Water Development Board 
is used with substantial moclif icatiens . 


1 *3 Organ! zati'on of the Report 

In Chapter 2 of this thesis, the development of the 
techniques of synthesis of streamflow is summarised and the 
description of the streamflow gonora-tion model adopted is 
given. The deecri.ption of the out -of -kilter algorithm is giver- 
in Chapter 3 along xvfth an illustrative example and the model 
used in the present study is presented. Chapter 4 contains the 
descriptions of the basins selected for study and the analysis 
of the river basin system. The results of the study a.re 
presented and discussed in Chapter 5. 
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1 • 4 Units Used in the Stiady 

In order to keep the study as close to a real prohlem 
as possible, actual data from published literature for certain 
river basins in India were used. These data are generally in 
P.I.S. system and hence the same units are retained in the study. 



OHAPTIH 2 


STHEMFLOW GSNSRATIOF 

2 . 1 Introduction 

A coronion constraint encountered in water resources 
management is the inadequacy of str'eamflow records. If policy 
decisions on design and operation of vrater resources system 
are based on inadequate information, the response of the system 
under a full range of conditions cannot be postulated. If the 
length of streamflow record is short, critical sequences of years 
of low and high runoff ' inherent in the statistical population 
of river flows may be missing. But no matter how poorly a 
brief record may identify the time frequency of years or seasons 
of unusually high or low flows, unless it is very short indeed, 
it will permit fairly reliable estimates of mean annual and 
seasonal flows and their variances. These statistical parameters, 
together with a few assumptions about the population of flows, 
can make it possible 'to construct statistical models that can 
generate synthetic flews of any desired lGngth(Maass et al. , 

1962 ). 

2 . 2 Ti me Series Components of Stream f low D ate^ 

From a statistical point of view, streamflow data can be 
regarded as consisting of four components (Kottegoda, 1970) 
viz., trend T.j., periodic or seasonal S^, correlation and 
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random compcncnts which can he combined simply as 

follows 


-ii-j. = 3^ + (2,1 ) 

A seqaonco: of values arranged in order of their 
occurrence is called a time scries. A time series is 
considered to be sta-tionary if the statistical properties 
characterising it are tiiae invariant. The non-stationary 
data can be made stationary by a simple transf ormation 
(McMohan and Mein, 1978). 

One characteristic of a time-series is persistence which 
relates to the sequencing of the da-ta.- In streamflow, per- 
sistence arises from natural catchment storage effects which 
tend to dGj.ay rxin-off; over a short period of time, high flows 
in one interval will tend to be followed bj^ high flows in the 
following interval. The longc'c the time period, the lesser 
the effect and for manjr streams it is negligible for annual 
f].ows (hcf-Iohan and Mein, 1978). 

The usual quantitative measure of persistence is serial 
correlation. Serial correlation coefficients may be calcu- 
lated for the correlation be. tween the flow in any given time 
period (for example, month or 3 '-car) and the flow in k time 
periods earlier where k is called the lag (k = 1,2, ...). In 
many studies, only the lag one serial cor’^-elation is con- 
sidered, that is, the persistence between an event and the 
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immedi? toly pro coding event. Lag one models have been 
shown to be operationally satisfactory in several studies 
(lotto ga da , 1 97C' ) . 

2 • 3 Jovelopment of Synthetic Btreamflow Generation Techniques 

itllen Ha son (1914) is considered to bo the first to 
recognise the desirability of extending hydrological data. 

He synthesised a rinioff sequence of 300 years by combining 
annual--mo3;n-~ilow series for 14 streams in which the floxfs of 
each were expr£;ssed in torms of individual, mean flows. This 
method in effect combines samples from different populations and 
is thus not preciselje applicable to any particular stream. 
Charles 3. ludler (1927) employed a deck of 50 cards, on each 
of which was printed a representative annus.l streamflow. By 
deailing this deck 20 times, he obtained an artificial record 
of 1000 years. The J^dsquacjr of this method depends on 
ho'w the vaJ.uea printed on the deck are determined. Further- 
more, the method has tne unrealistic limitation that the - 
largest flow in 50 years is also the lai-gcst fli.ow for the 
entire record. In this merhod, all periods of 50 years have the 
sane mean, the same standard deviation and the same range which 
is a major defect. F.B, Barne-s (1954) used a simil.ar method 
to that of Sudlcr except that the synthetic flows were appro- 
ximate}.y made norna3. variates xirith the same mean and standard 
deviation as the flows of the historical, record. Barnes 
introduced the improvement of using a table of random numbers 
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in synthesizing a 1000 year sequence of streamf loxv-s . The 
use of BariiGs method is limited to the representative annual 
flows of a single stream that arc approximately normally 
distributc'd and that do not exhibit serial correlation. 

Markov introduced the concept of a process in which the 
probebi3-ity distribution of the outcome of any trial depends 
only on the outcome of the directly preceding triaJ. and is 
independent of the previous history of the process. In 
this 'trial' is the passage of one year and its 'outcome' is 
the strearxiflow for that 3 ?'oar. If the probability distribu- 
tion of annua?.- s'creamf 3-Ow is either independent of previous 
atreamflows or corre3.ated with onl 3 '' the previous year flow, 
we have a 'simple' or 'lag one' ilarkov process. The concept 
has been extended to include cases of lag greater 
than one and the process has been tne basis of study and 
deve.lopments of streamflow generation pj^ocedures during 
early sixties (Maass, 1962). 

Brittan proposed the fo3.1owing Markov model to 
repi’cscnt actual streamf lows (iicr!oha„n and I-iein, 1978).' 

^'"i+1 = (1 - Y^)’‘ (2.2) 

where 

X. ,X. , = - annual run offs for ith and (i+l)th year, 

p = mean historical annual fl.ows, 
s = standard deviation of annual f3,ows, 

Y.. = annual lag one serial- coi'relation coefficient, 
' and 

t,. = normal random variate with mean of zero a.nd a 
variance of unity. 
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Tnxs mode], consists of two components % a deterministic or 
correlation component [li-T Y ^ a. random component 

f t « ( 1 — 1 

L \ 1 i ^ J j * 


In the anmial I-Iarkor model orecl.ined above, only two of 
tne four covoponents assumed to make up the stroamf lox^r process, 
as defined in equation (2.1) are a,ccounted for explicitly. 
Trend and periodicity are not considered. 

The most common form of periodicity relates to seasona- 
Ixty, particularly with respect to monthly flow generation. 
Here, the most .appropriate practical model is the one proposed 


Thome 


".na rxerxi 


( 1962 ). 


The algorithm for the Thomas 


Pioring season.al model i 


as 


-lows 




= pM ’■ 




(2.3) 


wnere 


P. , , ,X. 

XtI ’ X 


gunorated flows during (i+l)th and ith seasons 
reckoned from tne start of the synthesized 




D . 

3 


mean f3.ows during (i-j-l)th and jth seasons within_ 
a, repetitive annual cycle of seasons (if months 
p-rc being used, 1 £; j 12) 

least square regression coefficient for estima.- 
ting (j+l)th f.low from jth flow 


= Y. 


la±l 


(2.4) 



■I 


random deviate with 


1 1 



"o 


normal 


mean of aero, and 


varic.iice unity 

= staniard deviations of flows during {j+l)th and jth 
sessonvS, and 


= correlation coefficient between flox^s in jth and 
( 3 +l)th seasons. 


To use the model to generate nonthlj'' flows at a site, 36 
parameters - monthly means, standard deviations and lag one 
serial correlation ~ ai'e required. These are obtained from 
analysis of monthl,y historical flows. 


This modoj. is restricted to normally distributed flows, 
that is, t^ is considered to be a, normal, random dsvia.te. In- 
order to cater for non-noriTia.?- stresnf 3-Ows, the model can be 
.modified by any of the following alternatives. 


( 1 ) modify t., by an appropriate transformation 
(2) modify tnc strcamflow parameters and the model al-gorithms 
such that the final gor-e rated lata arc distributed like 
the historical flo-f upon which they are based 


(3) generate normal.ly distributed f3-0ws and apply inverse 
normr>3.izing equation.s. 

Ma.ta3.8s (1967) presented moment transformation equations 
vriiich theo?oGticall 3 r preserve the .moments and la,g one serial, 
correl.a.tion coefficients. This method assumes that the 
l.o.ga.rj.thiTis of the f3_ovrs aro norina.3.1y distributod. The procedure 



12 


ia firat to Galcul.ete e. series of logaritiims using a normal 
model and then obtain absolute flo-sirs by exponentiation. 

2 • 4 Method Used in the Present Study 

In the present study*, a method devolopod by Beard ( 1 972 ) of 
the Hydraulic Engineering Centre, U.3. Army Core of 
Engineers, for multi-sites and multi-periods is used. The 
following equations are given by him. 


X. 

1 .m 




A. 

1 


Si 


Si 


= V X. /F 
^ 1 ,m' 

m=1 


... ^ P 

2 (X^ ^ - X^ ) / (lT-1 ) 


m=1 


1 ,m 


= N 


S, ffl.m -Zi)h[(H-1)U-2)s5] 

m=1 ’ 


where 


(2.5) 

( 2 . 6 ) 

(2.7) 

( 2 . 8 ) 


X = logarithm of incremental monthly flow, 

Q = monthly recorded strcamflow, 

q = small increment of flow used to prevent infinite 
negative values of X for months of zero flow, 

X = mean logarithm of incremented monthly flows, 

U = total years of record, 

S = unbiased estimate of population standard 
deviation, 

g = unbiased estimate of population skew coefficient, 
i = month nimber, and 
m = year number. 
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I'acii inaix^idual f?LOv.' is then conTerted to ^ normalized 
stanaard variate, using the following approximation of the 
Pearson Type III distribution. 


t . 

1 ,m 




(2.9) 


Si 


^ ^ ^i.m /2)+1) 


1/3 


- 1 ] 


4 - 


Si 

r~ 


( 2 . 10 ) 


where 

t = Pearson Type III standard deviate, and 
iC = normal standard deviate. 


The above equations are used for generation. The generated 
normal standard deviates are converted back to flows by use 
of the following equations. 


'^i,m = - g./6)H-1) ~1 Wg. (2.11) 


'\r 

- 


= X- + 3. 

i,m X 


and 


where 


Q. = -q. + Mtilog Z. 


Q > ° 

i,in 


( 2 . 12 ) 

(2.13) 

(2.14) 


k' = monthly flow logarithm expressed as a normal 
standard deviate. 


2 . 5 Streamf3.ow Generati on 

Generation of streamflows is accomplished by starting 
with average va3.ues for all stations in the first month and 
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discarding the first two years of generated flows to ensure 
a rcal.ly random start. Maximuia, minimum f=nd average flows 
are obtained for the entire period of flows as recorded and 
for specified periods of years. 

Because of limitations in computer memory size and also 
due to increasing change of computational instability with 
larger matrices, the number of stations usable simultaneously 
in this raodel has been limited to eight. By including a 
few important stations from one set to the next set of 
stations, again limiting the total to eight, simultaneous 
flows for all the stations in a basin ia generated, preserving 
the important correlations. 



CHAPTBR 3 


0PTIMIZA2IQII MODEL 


3 . 1 Oeneral 

'^/’a.ter resources proolems are beconiing more complex 
and larger in size. Liieir efficient planning and design 
requires tiic.t the most powerful analytical techniques 
avail-sble be used. Simulation and optimization techniques 
coupled with high speed digital computers can, if used 
properly, pro-ride the planner -with -ra?-uable decision-aiding 
information. In this chapter, a brief description of the 
model used in the present stud.]’’ is given. 

''•2 lechniaues f or Outi miza tion of a Water Resources System 

Linear programming and dyn?imic prograimEing are two of the 


optimisation techniques that have been widely used in water 
resources system optimization studies. Linear programming 
has been used to solve many water resources problems by 
researchers like Rogers, (1959), Mar^lin ( 1962), Maass (1962)^ 
Tiioraas and itevelle ( 1956) and Loucks (1969). The use of 
dynamic progranraing for sol-eing wa.tor resources problems was 
pioneered by Hsd.I. and his co -workers (Hall and Buras, 1961; 
Hall, 1964; Hall and Howell, 1963; Hall and Roefs, 1966; Hall,, 
Butcher and Eso,gbue, 1368). Later, Hobasheri and Harboe 
( 1 370) and hut Cher and Sunder (1975) also used the dynamic 
pro,grajiuning technique for optiriiiza.tion studies. 
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In the prese-at 
makes use of the out 
Fulkerson (l96l) is 
have demonstrated th 


study, an optimization technique which 
-'■of kilter algorithm (OKil) introduced "by 
used. Iiiimael.blau end G'laoghaire (l974) 
e use of this algorithm in solving a 


water resources optiraization prohlem. 


7-3 Ou t-of-hilter AlA'Orith m 

This algorithm makes use of the concept that the vrater 
resources system can he represented bj,!' a series of nodes and 
arcs in a capacitated network form analogous to electric 
circuit. The prob?.em solved 'bj OKA is essentially a linear 
programming problem with the special feature tha.t a number of 
equality constraints exist. Computational results for some 
large scale problems show OKI to produ.ce a solution in one- 
twentieth to one fiftieth the tirao of standard linear progra- 
mming codes (Himmelb3.ou and O'Laoghaire, 1974). This is 
due to the following reasons t 

1) all oporations are additive (i.e., no multiplication or 

decision takes place) 

2) no matrix inversion is necessary. 

The use of OKA to solve a water resources problem is 
illustrated by an example taken from Himmeiblau (1974). The 
probl-cm sel-ected is a minimum cost circulation problem. The 
various symbols used in this example are shown in Tab3.e 3.1. 
The iiiiie possible mutually exclusive 'kilter conditions' for 
each arc as the algorithm proceeds to seek an optimal 
solution is shown in Table 3.2, 
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Table j.] Symbols Ised in iiie Out-of "■Kil.ter 
llgorithm 



-O 

± 


ij 





= Benefit of pa.ssing one unit' of flow through arc (i,j) 

= Ploxi in arc (i, 3 ) 

= Lower flow capacity of arc (i, 3 ) 

= Total cost of transporting one unit of flow from 
node i to node j 

= Upper flow capacity of arc (i^j) 

= Price of one unit of water at the node 3 ’ 



Tab le 

3.2 

Possible kilter-conditions 

for an arc 

State 


q- -! 


In Kilter ? 



3.3 



A 

1 

V 

0 

f = 1 

Yes 

B 

1 

= 0 

1 f U 

Yes 

C 


< 0 

f = U 

Yes 

A. 

q 

> ° 

f C 1 

ITo 

h 

q 

= 0 

f < 1 

No 


q 

< 0 

f ^ U 

No 

Ag 

q 

> 0 

f > 1 

No 

Bp 

q 

= 0 

f > u 

No 

C\1 

0 

q 

A 

0 

f J> tl 

No 
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Fig. 3.1 shoifs the network representa.tion of the example. 

Tne ordered triple (1^^, -; ) is shox/n on each arc. The 

original arc f.loxf ana node numbers (% values) appear as 
underlined numbers. Arc (1,3) is out of kilter and in state 
Ap. This arc can be brought into kilter by changing flows in 
the closed path (1,2), (3,2), (l,3), by increasing flows in 
forward arcs and decreasing flows in reverse arcs. The process 
is continued till all the arcs are brought into kilter. The 
node numbers (re values) keep changing with each iteration and 
a labelling procedure is to be followed. 


3.3.1 labelling procedure 


1 . 


2 . 


If an arc (i,j) in state A^,B^ or is brought into 
kilter by increasing the floxfs, then the node j is 
labelled (i"*”, e(j)). This means that node j may receive 
e(j) miits from node i. e(j) is taken to be 
the arc is in state A.^and ( ’'‘'•j j_ j ) state or . 

If an arc (i,3) in state A^jB^ or Cp- is brought into 
kilter by deex-easing the flows, then the node is labelled 
(3”, e(i)) meaning that the flow from node i to node 3 
can be reduced by e(i). e(i) is taken to be (f^^-lj^^)if the 
arc is in state Ap and (f^^-u^^) if in stafe B2 or Cp. 

An arc in any of the possible remaining states A,B or 


C is in kilter and its flow should not be changed. 



INDEX 




(a) Example problem (b) Flow values and node numbers in iterations 2,5 & 7. 
Fig. 3-1 Minimum cost circulation problem 
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3 . 3.2 So lutio n 

Suppose tnat arc (i,j) is out of kilter and node i has 
teen labelled. Now, the aim is to find a flow augumenting path 

from i to j in such a way tha-t in-ki 3 .ter arcs in the path will 

« 

not be driven out of kilter. An arc that originates or termi- 
nates at a la.belled node is considered and an attempt is made to 
label the node at the arc's connecting end. If node j is 
labelled, a flow augumenting path has been found and the flow 
in connecting cycle is changed according to the label value. 
Thus arc (i,j) is brought either into kilter or made less 
out of kilter. Now, another out-of-kilter arc is selected 
and the procedure is repeated. 

The five principle steps of the OKI are given below. 

1) Find an out-of-of kilter arc ( 1 , 3 ). If none, the optimal 
solution ha.s been found. 

2 ) Determine whether the f3.ow in tlio arc should be decreased 
or increased to bring the arc into kilter. 

3) If the flow in the ai-c is bo be decreased, find a path 
from i to 3 along which "the flovr can be increased without 
causing any arc to become out-of -kilter. Increase the 
flow in the path and decrease the flow in ( 1 , 3 ). If 
(i, 3 ) is now in kilter, go to step 1 . If (i, 3 ) is out of 
kil.ter repeat step 3. If no path is found go to step 5- 
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4) If ‘cne flovr in tiae path is to be increased, find a path 
from 3 to r along which the floxAr ca.ii be increased without 
causing a^ny arc to become out-of ~kilter . Increase the 
flow in the path and also in (i,j). If (i,j) is now in 
kilter, gc to step 1. If (i,o) is out of kilter, repeat 
step 4. If no path is found go to step 5. 

5) Change the % values and repeat 2 for arc (ijj), keeping' 
the same labels on all nodes already labelled. If the 
node numbers become infinite, no feasible solution is 
possible . 

Following these five steps, the ezample problem is solved. 
Seven iterations are necessary to obtain the optimal solution. 
The results of the first iteratioji is shown in Table 3.3. The 
final results, vmen all the arcs are brought into kilter is 
shOT-m in Table 3.4. The Fig. 5.1(b) siioxfs the results of 2nd, 
5th and 7th iterations. 

A small netx^ork x^ith 4 nodes and 7 arcs takes 7 iterations. 
In complex water resources system, where the number of nodes and 
arcs are larger and where steps 3 and 4 of this article are 
repeated, the number of iterations xm.ll be enormous and a high 
speed digital ccmputer becomes essential to solve such problems. 
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TaloD-e 3.3 Results of Iteration 1 


Arc (i,j) 


-TC - 

D 

-b. . 

13 


^13 

State 

In kilter? 

(1,2) 

0 

0 

2 

2 

0 = 1^2 

1 

Yes 

(1,5) 

0 

0 

5 

5 

2 

A 2 

Ho 

(2,3) 

G 

0 

1 

1 

0 = I 22 

A 

Yes 

(2,4) 

0 

0 

3 

3 

^ ^24 

Aj 

Ho 

(3,2) 

0 

G 

1 

1 

2 > 1^2 

Aj 

Ho 

(3,4) 

0 

G 

6 

6 

° < I34 

A, 

Ho 

(4,-1) 

0 

0 

0 

G 

2 < I41 


Ho 



Table 3 . 4 Re suits 


( i , D ) 

11 . 

X 



,2) 

0 

-4 

2 

,3) 

0 

-5 

5 

,3) 

4 

-5 

1 

,4) 

4 

-1 

3 

,2) 

5 

-4 

1 

,4) 

5 

.,.7 

6 

,1) 

7 

-.-4 

0 
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Iteration 7 


^ij State In kilter? 


•2 

C\J 

11 

<\I 

0 

Yes 

0 

1^2 - ^ - ^12 

B 

Yes 

0 

0 = I23 

B 

Yes 

0 

^24 - ^“^24 

B 

Yes 

2 

0 

11 

H 

ro 

A 

Yes 

4 

1 - 134 

A 

Yes 

7 

3 = 1^^ 

A 

Yes 
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'3.4 -Jescrip'tion oT t.he Model Jsed 

in Ain.o ontion Models Al- 3 ? devel.opecl "by Texas ¥ater 

development Board (T¥]jB-2i, 1972) is used. inis is a general 
hydrologic model of surface water resource sysee.ea. The model 
is designed to ansJ.yze the simulated muD-tiperiod operation of 
any interconnected configurations of reservoirs, pump canals 
and pipe lines on a steady state monthly or seasonal basis. 
Substantial modifications are made in this model to suit the 
present requirements and also to improve the efficiency of 
the model in redistributing the deficits evenly to various 
demand nodes. 

5-4.1 Math em atic a l model 

The out -of -kilter algorithm, OKA (rlimmelblau, 1974) used 
in the model solves the following problem i 

For the arcs ajia nodes defined in the system, 


minimize 

CO st 

" E E °lj 

1 i 

(3.5) 

subject to 

q. . 


(3-6) 


a . . 

-X3 

^ hi 

(3-7) 


Sq. . 

"ID 

- • • =0 for each i 

(3-8) 


and 


q . . 0 for each a. • 


(3.9) 
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where h . . = cost of passing one unit of flow through arc a.., 

J i J 

q . . = cpaantity of flow passing throiogh arc a. ^ 

J 1 J 

H. . = upper bound on arc a. 

L . . = lower bound on arc a. 


13 


13 


and 


a. . = a.rc connecting node i to node 3. 

J- J 


Eqn.(3.8) when applied to the various nodes in the network 
gives the following equations. The various terms used in the 
equations are defined in the Table 3.5. 

(1) Initiad. storage node ; 


n 

r 

S 

3=1 




(3.10) 


(2) Input node % 


m jd 

T r 
S S “jk 


k=1 3=1 


= 2^-5 (Inflow from Source Rode) (3.11) 


(3) Demand node ; 


T 


n 


n^ 


k=1 

( 4 ) Import node ; 

n. 




(3.12) 


(Outflow to Sink STode) 


T "t 

2 S ^3^k from Source Node) (3.13) 

k =:1 3=1 
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( 5 ) Spill node s 


2 2 ^3k 

k=1 3=1 


. (Outflow to 
Sink Hode) 


(3.14) 


(6) Pinal storage node 


n 


S 

3=1 




4:4; 


= X„ (Outflow to (3.15) 

Sink Node) 


(7) Net Dalance node 


nr , :r 

-'"i 


iv. ^ 

a m 


= 0 


( 3 . 16 ) 


(8) Reservoir nodes t for the existing arcs in the net\'70rk. 


N 

2 

i=1 


2 ^i3k ■■ 


N 

2 

i=1 


2 ‘^jik 


O' . P . I 

3 3I5; 


nt nt 


+ “3k ° 


(3.17) 


3 = 1,2, . . . , n^ and k = 1,2, . . , , T . 


(9) link junction nodes ; for the existing arcs in the network. 

¥ 


N 


Q- ■” v. Qmv “■ + 6.-li_ + - D_.i_ = 0 


2 '“^i-jk " 2 '=^3ik - - 3 " 3k ^ "3"k '■ " 3 k “ "jk 

i=1 i=1 (3.18) 


3 = n. 


r-^-l * ■■‘'r+2 


. . . , n^ and k = 1,2, ... T, 
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Table 3.5 

Dei iiiitioii of Terms used in AL ~3 Model 


Network flou pro b lo rn 

0 

b z= Coot' 01 moving one- unit ox floxAr from node i 
to node j 

1 . . = lower bound on flow from node i to node j 

-L. J 

N = Number of nodes in the network 
q. . = Plow from node i to node 3 

j- j 

U. . = Upiper bound on flovr from node i to node 3 

1 J 


Node ba 3 .anc; 


lations 




■k 


= Demand at node 3 in time period k 

c 

= Rate of import of water in time period k 



= Rate of spill from reservoir or link junction 
3 in time period k 



Plow between reservoir or link junction i and 
reservoir or link junction j in time period k 


S.-, = Storage contents of reservoir j in time period k 

3 k 

a ., = Input to reservoir j (unregulated inflow) in 

^ time period k 



= (Demand + evaporation) from reservoir j in time 
period k 


6 . - 1 , if j is an import node- 

3 

=0, if j is not an import node 

e . = 1 , if j is a spill node 

3 

=0, if j is not a spill node 
At = Unit of each time period 


units®' 

X(l^/t) 

1^/t 

1^/t 

1^/t 

iVt 

1^/t 

1^/t 

l^t 

1^ 

1^/t 

1^/t 


contd . . 
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Table 3.5 (ccntd...) .. 

Subscript s 

i , j = Nodes 

Ir. = Time period 

= Number of reservoirs 

= Npunber of reservoirs plus link junctions in each period 
T = Ntiniber of time periods 

1 and t are used to designate xmit of length and time respectively 

(1 /t = volume per unit time) 
b 

X is given in relative cost units 

c *■ 

only ozie import node is considered. 


mv 


he above model' is used for the sj'-stem of Basin BA. A des- 


cription of this basin is given in the next cha; 


p b • 



OHIPTSH 4 


ilTALISIS OIr reiiE MODEL 

4 • 1 I ntro due ticn 

Two major river basins in Indie are te.ken for study. 

Both the Streamflow Genera,tion model and the Allocation model 
are applied to one of the basins BA and only Streamflow Genera- 
tion model to the other basin BB. Basin BA has inflow data 
varying from 8 to 38 years at various nodes while for the basin 
BB , inflow data, are available for an almost constant period of 
ton years. The hj^drology of both these basins are governed by 
monsoons. 

4 • 2 Description and Modelling of Basin BA. 

The system consists of 15 reservoirs and 4 diversions and 
they are connected by river reaches. A description of this 
basin is given by Rramamurthy (i960). The same basin is analysed 
here with 3..onger and more reliable data. 

The rj.ver system of Basin BA in the network form is 
shown in fig. 4.1. Bodes represent reservoirs and diversions 
while tne links represent the ri.ver or canal reaches. The 
details of the reservoirs and diversions are shown in Table 
4.1. The details of power stations are shown in Table 4.2. 
Table 4.3 gives the details of system link connections, the 
bounds of the links aaid relative costs. The upper bounds ox 



1 

<jlO| 

t 


2A 


< 




i 

A 

O 



Power Generation 
Reservoirs 
Link Junction 
Artificial Node 
River Reach 

Power Canal 


Fig. 4‘1 Spatial network of basin BA 
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Table 4. 1 L’Gtails of reservoirs and non storage 

nodes for Basin BA (values in 10 M.Cft) 


Node 

No. 

Reservoir 

Maximum 

Capaci tv 
Minimum 

Annual 

demands 

Mean 
f 2.0W 

75 depen- 

dable flO'?>r 



Stora.ge nod 

es 



1 

30800 

15800 

2532 

3397 

2741 

2 

3370 

560 

4240 

8477 

5947 

5 

3765 

717 

14840 

2397 

2104 

4 

51 16 

31 1 

8835 

9374 

7319 

5 

4224 

1182 

7753 

30909 

27121 

6 

3774 

713 

4094 

4799 

3744 

7 

2820U 

1 9400 

28086 

2680 

2194 

8 

11136 

6000 

8125 



9 

955 

1 80 

7038 

16244 

12777 

10 

7154 

850 


10132 

7883 

11 

13300 

1560 

21760 

1887 5 

14798 

12 

9314 

423 

- 

13116 

11277 

13 

3960 

940 

600 

8305 

6556 

14 

1 1264 

781 

8531 

C\J 

O 

22931 

15 

3640 

280 

654 

12093 

9113 



I'fon~stora,g. 

i nodes 



16 , 



11279 

10359 

6990 

17 



30427 

36421 

28380 

18 



4200 

11494 

10085 

19 



1150 

22750 

16292 

20 

•m.*. 

~ Artif icia.l spill 

-- 

- 




node 





Artificial 

power nodes 



21 



5990 

- 

- 

22 



4300 


- 

23 


-H. 

6750 


— 

24 



■2560 

- 

— 

25 



4259 



26 



4930 



Bor 

tte.- full 'basin 


192933 

252754 

223007 



Table 4.2 Details of power stations in Basin BA 
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rP fP 


P 

cd 


cd P^ 
o 


CO 

0 

p 

d 

H 

o 

d 

H 


Cd \n o p 


Power station bctuivcalent to all u/ s power sta,tion 
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Table 4 -.:^ l)eto.x!Ls of link connections for 
Basin BA 


Link 

No. 

Node 

Prom 

To 

Maximum f lOi^r 
in to M.Cft 

Minimum flow 
10 M.Cft 

Relative 

1 

12 

5 

50 , 000 

0 

5 

2 

2 

5 

50 , 000 

0 

5 

3 

13 

5 

50,000 

0 

5 

4 

4 

5 

50,000 

0 

5 

5 

22 

16 

710 

0 

2 

6 

5 

3 

50,000 

0 

5 

7 

6 

3 

50,000 

0 

5 

8 

3 

1 

50,000 

0 

5 

9 

1 4 

8 

50,000 

0 

5 

10 

8 

9 

50,000 

0 

5 

1 1 

26 

Q 

50,000 

0 

5 

12 

1 1 

16 

50,000 

0 

6 

13 

7 

20 

99,000 

0 

5 

14 

10 

1 1 

50,000 

0 

5 

15 

9 

1 

50,000 

0 

5 

16 

16 

1 

50,000 ■ 

0 

5 

17 

1 

7 

99,900 

0 

5 

18 

18 

5 

50,000 

0 

10 

19 

19 

1 1 

50,000 

0 

10 

20 

20 

17 

99,000 

0 

4 

21 

15 

26 

50 , 000 

0 

5 

22 

1 , 

■1 

7,775 

0 

0 

23 

12 

■ 5 

550 

0 

1 

24 

10 

21 

790 

0 

0 

25 

1 1 

22 

710 

0 

0 

26 

12 

23 

722 

0 

0 

27 

13 

24 

320 

0 

0 

28 

14 

25 

470 

0 

0 

29 

15 

26 

590 

0 

2 

30 

10 

1 1 

830 

0 

I 

31 

11 

16 

530 

0 

3 

32 

7 

20 

9520 

0 

0 


34 


ca.nal lixiks are equal to the inaximuui quantity of ua.ter that 
can ht' passed xhrougn turbines of power stations concerned* 

4.2.1 Demands 


111 India, irrigation projects are to be designed so as to 
meet /5 percent dependability and power , pro jects to meet 
90 percent dependability. In the present study irrigation demands 
have priority compared to power demands. The annual 75 percent 
dependabl.e flows are reallocated to the various nodes in the 
system excluding the evapora,tion losses where evaporation losses 
are known. There are a number of minor schemes in the basin 
and their demands are added to the appropriate nodes. When 
the inflows during a year are less than, 90 percent of the 
75 percent dependable inflows, that year is termed as a dry 
year. I'or such years, demands are reduced to 90 percent of 
the planned dumands. Similarity, when the inflows are greater 
"by 10 percent of the 75 percent dependable flows in a year, 
that year is considered as wet year and demands are proportiona- 
tely increased. The demands for the monsoon months are kept | 
at average level since no a.:;ditional reservoir storage space is 
available during these periods. The average year demand for all i 
nodes is snown in Table 4.4. The monthly evaporation rates for 
the reservoir are shown in Taole 4-5. | 

The monthly reservoir water spread areas and the heads I 

acting over the turbines for the throe hydrological states j 

is sh'-.'X'jn in Tcodcile 4.6 a.nd Ta^ble 4.7 respectively, : 



Table 4.4 Average Year Demands at Various lodes for Basin BA 
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Table 4.7 Average Monthly Povrer Head Rules 
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4.2.2 Determinat ion of aub b.q..q_ir, yi eLas 


'Tne basin is subdivided into tvrelvG sub basins. From the 
strearaf lovf data and cletai3.s of withdrauaJ.s available in each 


o:t tho sub basin, the raean yields of each sub bg,sin is 
e stijiia tod . iin illustrative example is given below on the 
method advOpted to obtain the mean yield of one of the sub 
basins. Tho sub basin chosen is shown in fig. 4.2. 


Catchment area of the sub basin ; 6939 sq. miles 
Infl.ow details available i 

SI. Stage discharge Catchment Mean 


No. 


site 

sq. miles 

yield TMC 

1 


SDl 

2134 

173.40 

2 


SD2 

761 

112.16 

3 


s:d3 

907 

89.19 


Estimation of sub basin yield ; 


Yield for 2l34 sq. miles upto SDl 

r= 

173.40 

TMC 

yield for 761 sq. miles upto SD2 

= 

1 12.16 

PMC 

Yield for 907 sq. miles upto SD3 

= 

89.19 

TMC 

Estimate yield for T4 basin for^ x 
907 sq. miles as per SD3 (both 
in the same rainfall zone) y 


89.19 

TMC 

Diversion at Pi 

= 

86.98 

T!4C 

Withdrawals at P2,P3 and P4 

= 

2.60 

TMC 

Estimated yield for remaining areax 
of 2230 sq. miles at the rate of J 
18 M.Cft/sq. miles 5 

=: 

40.14 

TMC 

Mean yielcl of the sub basin 

say 

,590, 

593.66 

Tie 

TMC 




SDi ,5D2*--- stage discharge sites 

ft . 

Fig. 4-2 The sub basin SB A-1 



liic oc.iiAc procGaiLCkj is aaopteci to ca,lc'iilatG tha mean 
yield of al], the twelve sub basins. The mean yields of 
various sub basins are shown in Table 4.8. 

4-2.f) S tr c anif low Igrjjrp -h nn 


4 . 2 . 3 . 1 I)at cO B a so 


Streamf lows of sufficiently long periods are needed to 
siniula.te tho system, for large river ba^sins, it is difficult 
to collect the relevant concurrent streamflow data needed at 
representative points. Reliable measurement of data during 
the periods of peak flow is difficult and expensive. For the 
ba.sin BA, gauged data are available for periods varying from 
8 to 38 years, the longest being available a,t node 4. 

li’or the present study, the virgin flows (the flow 
reaching a. node from the independent catchments between the 
upstream nodes and tho present node) arc required at all the 
nodes of the system. The virgin flows are estimated at the 
various nodes by the follovring methods. 

1 ) The procecLui’e adopted to determine the mean annual flows 
for the sub basins as explained in a,rticle 4.2.2 is 
adopted to estimate the mean annual inflow's at all nodes. 

2) These mean flews are proportionately altered to maintain 
the total basin mean flovr to be equal to the total mean 
yield of the whole basin. 
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Table 

4.8 

sub 

Yields of 
basir^of Basin BA 

various 

SI. Ho. 

Sub 

basin 

Mean yield 
(T.M.C.) 

1 

SBA 

~ 1 

590.00 

2 

SBA 

~ 2 

60.00 

3 

SBA 

- 3 

180.00 

4 

SBA 

- 4 

70.00 

5 

SBA 

” 5 

460.00 

6 

SBA 

- 6 

70.00 

7 

SBA 

~ 7 

230.00 

8 

SBA 

-8 

540.00 

9 

SBA 

- 9 

85.00 

10 

SBA 

- 10 

65.00 

1 1 

SBA 

- 11 

18.00 

12 

SBA. 

- 12 

80.00 


48 


5) 


lor '..".16 i'iod.6 'Ui'-Ci6i' co!n.si.d.6!r2,'t;.1.on 
gauging station whicii i,s close-st 
to represent the dispersion of f 
mean monthly flows. 


, the fT.ows at the 
to the node are .assumed 
lows 67i-Gh. re-spect to the 


4) The annual mean flow at the gauging stations nearest to 
the node is proportionately altered such tha.t this mean 
is equal to the mean calculated in step 2. 

5) 'The effect of upstream releases, if any, are neglected 

in this estimate. The variation in the dispersion caused 
in the flows due to these upstream releases is also 
neglected. This simplification does not significantly 
affect t.hc results of simulation since, the total meanflow 
at t.he node is maintained. 


Ta.ble 4.1 shows the annual mean flows for the various 
nodes. In the recoiistitutod flov/s at various nodes, nega.tive 


values are asooigned if the flow at any month is missing. 

The so are .rec'" gnized by the streamflow generation model and 
proper values are e-stimated to these points by multiple re- 
grooGion. This forms the data base for streamflow generation. 

4 . 2 . 3 . 2 Strea mflow Generation of Baain .fBA 

There are 26 nodes in the model of which nodes 20 to 26 
are artificial nodes and hence virgin flows are all zero at 

I . . 

those .nodes. Due to limitation cf computer memory size, only 
8 nodes can be handled at a time for reconstituting the missing 
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fl- .Vo 1 a ting xhc streanf lows. Hence; gr-oups of 

n..cLo& C..O.,, fc:',r._fcn. at ..a time. 3o;oe impcootant nodes are kept 
cohmcn f. jo (?ocii sot of generatiovi of flows to aeintain the 
cro.os c :aTo3,.ation aio^ng various no4es. The n^:.ioo are grouped 
into four sets. The first set ccnsists cf nodes 1 o;., 1 2 ^ 4 , tC , 
9,16,1 £inci / . These a,re representative ncdos taken from the 
upstroaiti, ciidnle 8,nd downstreaiui pcrtions of the river hs-sin. 


Prof.'i those, no dCvS 4,9,12 a.nd 14 are carried fcrw8.rd and 
node)s 2,3,15 end 18 are added. For the third set, nodes 5,6 
a, lid 11 aro cciisidered a.lcng with 4,9,10 and 16 from previous 
sets, .e'-^r the last set, nodes 3,7 and 19 are chosen with 
n„dGs 1,6,9 and 10 from the earlier passes. This forms the 
schemG f-'r generation of synthetic flows. 

Streamf Ico'fs for 138 years were generated f-or all the 
nodes including the reconstitution of missing data in the 
38 years of base period. 


4 . 2 . 3 . 3 Compa r ison of generated flows 

The goner a ted streamflcws are statistically analyzed to 
deternino their, mean and standard deviations. These are 
shown in Table 4.9. The values cf the statistical parameters 
for bv th rc- corded a.nd reC'-instituted flows as well as for the 
generated flows compare favourably. The statistical parameters 
for the ouitire basin is shown in Table 4.10. 
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4 . 2 . 5 . 4 3i nulat ion ^ the system 


J.hb; llJ.o cation mcdel is used to siuulets the system. 

To reduc:.,; Lht. me.norj/- size in the computer, a year is divided 
into live iiOii iji'ily periods and a last period consisting of 
the 1 t;i,!.ai.aing montxis of the year. The nanimun size of the 
network go he used consists of four years. Initia.l reservoir 
contents should be l8,rge enough to see that the simulation 
starts) with a feasible solution. 105 years of generated flows 
arc used for simulation. 


Various priority values, benefits for the reservoir 
storage arcs and costs for the river and link arcs to 
achi'iVe proper operation of the system are decided on the basis 
of y, few pK’oliiainsiry computer runs. The simulation is done 
using 33 years of recorded and reconstituted flows so as to 
reduce tlu; cost of computation tine for the basin and the 
perfornance of the system studied. The final costs used for 
unit in tho river and canal arcs arc shown in Table 4.3. 

The rv,: cults arc presented and discussed in Chapter 5. 

4 * 3 Descri pt ion and Modelling of Basin BB 

Tho river system modelled in the network form is shown in 
fig-. 4.3. All tliG reservoirs and diversions are shown as nodes 
and river reaches as links. The system consists of 21 reser- 
voirs and 6 link junctions. Reservoirs 4>5jlO,12 and 14 
have poW'-i’ sta/tion below the respective dams. Provision is 





(26/ A Reservoirs 

X O Link Junction 

j }'j Power Generotipn 

— River Reach 

/ Power Canal 
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made st tiiese reserYOirs to pass the water thrcugn the power 
stations il tiio irrigation demand ig low and this water flows 
to the aowne;cro.Si-:i. reservoirs aftir power is generated. Reser- 


voir 9 io an equivalent reservoir roprosentiin 


all the up- 


,sbrcd.ii u. fci li .-.c',tio.d ot water. j?roi;i r_‘.aervoir 1 water is 
diverted to two reservoirs viz. 10 and 20. Similarly, water 
is diverted to reservoirs 15 and 17 from reservoir 16. 


Rci: orring to the fig. 4.3, nodes 1 to 22 form the system 
nodes. iodoB 23 to 27 are added for calculation of power 
dovelojied at various power stations. 


Binco only streamflcw generation model is applied to 
thlE? oystivi, d: tails other than that required for this model 
are not .■';i.vuri there. 


Tht:; vn'im yields of the various sub basins of this 
basin wstiiuati,'d on the same way as that of basin B1 are 
sh.o''Wii in To bit..' 4.11- 


4 • 3 . 1 Streanilow gwneraticn of Basi n B3 

Thu procedure adopted is same as that followed for the 
bacin ,BA uxcopt that in this case 5 ■■^ots of stations are used, 
Thu first out consists of nodes 1 ,5,354,12, 11 , 14 and 21. 

Trviu thusu riCvbns 1,5,3 and 11 are carried forward and nodes 
2,7,10 and 16 arc added to form the second set. The third 
■set C:'.jnsi.;-;it 3 ''''f nodes 5,12,7 and 10 frora the earliu-r sets 
to,:.,th„r «ltb nobos 8,9,20 and 6. Ifodes 18,15 and 13 together 
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Table 


1 1 


of Basin BB 


Yields of sub ba.sins 


Sl.iTo. 

Sub 

basin 

Mean yie 
(T.M.C 

1 

SBB 

~ 1 

179.49 

2 

SBB 

~ 2 

76.44 

3 

13BB 

~ 3 

101.84 

4 

SBB 

~ 4 

171.39 

5 

SI;B 

- 5 

198.91 

6 

SBB 

^ 6 

35.36 

7 

S.HB 

- 7 

'175-59 

8 

SEB 

- 8 

266.60 

9 

b:i?b 

- 9 

812.83 

10 

n'VbX'j 

k)J 

- 10 

786.06 

1 1 


,-11 

706 . 92 

1 2 

g;3B 

12 

419.21 
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with nodes ^, 5,20 and 6 forn the fourth set. The last set 
consists ci nodes and 20 from earlier sets along with 

nodes 17,19 and 22. 

The length of data available varies fron 6 to 10 years 
and they a.rs all concurrent. 

The results of the flow generation are shovm in Table 4-1 
and it is found that the recorded and reconstituted data and 
the generated data agree statistically. The statistical 
parameteiB f or the whole basin are shOTca in Table 4.13. 



CD 

rH 

cd 

o 

•H 

u 

o 

-p 

CQ 

•H 

Ch 

O ^ 

O'!) -p 
0 Vl 
•H O 

-p . 

0 

P 40 

o ^ 


Jh 

rH 

cd 

o 

•r-j 

-P 

m 

•H 

+3 

+3 

m 


H 

0 

?H 

05 

m 

Is 

o 

rp 

pq 


OvJ 

0 

r-j 

.Q 

05 

ei 


p 

! 03 LH 

gc 

' §^5 

PQ 

«! 

•r 4 

i 

I 

m 

H --4 

c6 

-H trs! 

P 

kf — 1 



o I 

1 

<H 1 

1 

W 

^ <M 

Is 

^ 05 —1 

o 

H 




•x:} 


0 

rO -r- 

-P 

0 

0 

k 


0 


0 

# 


Oj 

a f 

CD --r 
1 -^ 


O Ch* 


O 00 


o H 

O' 


+ 3 ^ 

Ph'O 
0 
CO 


<tj. 


i H ' 

f 


0 


Is ^ 
P C\l| 


H 

cdi 

> 

q{ 

O) 

0 , 


n 


CD 


!> 00 
CM Lf\ 

o o 

IS^ 

CO m 
03 

cr^ o 

LOv 

8006 

8458 

15668 

17002 

1 1828 

12199 ' 

O KO 
EA in 

71 

71 

- -- 

182 

247 

-- o 

O 

VO ITi 

i>- r- 

90 

100 


120 

151 

5 

14 

58 

68 

o — 

1 10 

124 

MO 

00 o 

90 

107 

95 

108 

CM 

K\ tr\ 

93 

97 

10 

12 

108 

125 

CM tr'^ 
r>“ CM 

82 

87 

17 

17 

1 40 

146 

15 

17 

1 1 1 

120 

^r^ ^ 

(T) sr\ 

CM 

114 

128 

£> O 
CM iSO 

261 

285 

rrv ^ 

145 

160 

710 

955 

o m 

LO MO 

50 

57 

CM VD 

00 cn 

^ CT\ 

K\ K\ 

o ^ 

CM CM 

ro CM 

00 CM 
m MO 

iXN CO 

VD {>- 

76 

132 

1495 

1564 

VD CM 

m o 
m 

l^l 

199 

2810 

3216 

1694 

1840 

297 

508 

l>- m 

O 'xO 

CT) o 

1984 

2065 

171 1 

1867 

5992 

7220 

o o 

r>. 

CO O 

KN 

00 VD 

O VO 

O Ir\ 

4168 

4089 

1993 

1890 

^^^ VO 

OA 00 
ir\ v£> 

CM CM 

CO ^^^ 

O 

CM C- 
CM 

O 00 

^ V£> 

00 CO 
^o. |>r\ 

937 

1810 

1649 

1521 

1271 

1320 

1114 

1105 

1 604 

1425 

1 181 

1099 

O ISO 
!>- 00 

CIS 

o 

MD D 
'O vD 

o 

M- a. 

996 

950 

CM iX'. 

(T) 00 

CO 

in'' CM 

{?Ov 

107 

161 

W CtJ 

g cis 

W d 3 

p :1 d 3 . 

fjj ' c !:5 

IP ■ Ci 3 


6 


contd 


Table 4.12 (coatd. 


o 

00 


LA 

CD 

CA 

VO 

O 

I ^ 

CA 

CA 

{> 

1 



CM 

1 

1 

tA\ 

{An 



— VD CO 


CM { 


00 O 


m 


t O'! 


CO 


o 


VO 


LA 




CM t 


00 in 


O 

ro VO 


CM O 
LO O 
^ LTV 


ir\ (T\ 
O i>- 
LPv 'xh 


i>- tS^ 


L^^ VO 
VO Ln 
Od CM 


LA 'vO 
VO VO 


ej 


'M- 


Cn LA 
LA, VO 


CT\ ^ 
^ LA 


LA 

CA 


Aj ^ 

o ^ 

KA LA 


LA VC 
LA VO 
MO LA. 


Cb' 


Ch LA 
LA 00 

^A K'x 


O 


in 

-r— 

VO 

00 

LA 

CM 

CA 

LA 

VO 

!>- 

O'N 

EA 

00 

{>- 

' -o'" 

^A 

CA 

VO 

00 

(A 

LTn 

•M" 

{> 

CD’ 

At 

CM 

VO 

VO 



LA 

LA 

VO 

CO 

LA 


lA 



LA 




O 


CA O 


Cvl tA 


O CM 
tA LA 


LA ^ 


CM CO 
CA CA 


— 00 
VO VO 


CD 00 


CA 

a o 

CM -r- 


VO ITS 
O CA 

CM -- 


CA 
VO O 
LA EA 


C\J 

CAi 


CO 

G\ 


CA CM 
CM ^0 


VO CM 

C\i iA\ 


^ {>- 

C'l CA 


Lt"'\ 

C-'- 


VO 

O 

O 

A- 

CsJ 

!>- 


O 

o 

A 

Al 

CM 


LA 

CO 

CM 

LA 

CM 

CM 



EA 

LA 

A- 

VD 

LA 

EA 


VO LA 
LA ^A 



LA 

CM 

o 


VO 

CM 

00 


LA 

00 

LA 

— 

IA\ 

O 


o 

CM 

CM 

CO 

O 

>*— 

00 


C^i 

LA 

A? 

LA 

LA 




CD 

LA 

VO 




VO 

O 


LA 


LA 

O 

LA 

V— 


o 

o 

LA 

A 

LA 

LTn 

'M' 

Al 

CA 

o 

CA 

00 

o 

VO 

LA 

LA 


LA 



o 


A- 

VO 

LA 







A 

A 



O- sj- 
CM CM 
VO VO 


H-Ht r-h 

j,-L{ 


cb 


LA CM 
A- EA 
CA ■ O 
LA 


62 


o 

LTV 


VO 


00 {> 


CA CA 
'M" 


O tA 
LA LA 


VO VO 
VO VO 


■.J 

D 


o 

o 

LA 

VD 

'M" 

O 

LA 

CA 

CA 

VO 

-vb 

CA 

CA 

o 


'M" 

LA 

LA 

O 

O 

O 

vM- 

'T— 

.r- 


CA 

A 


VD 

CO 

CD 

LA 

-r— 



A 

'■M“ 

00 

£>- 

O 

O 

o 

VO 

LA 

(A 

LA 

LA 


o 

> 1 — 

A 

A 

EA 

LA 

VO 


LA 

VO 





A 

A 

'M" 

LA 

A 

A 


, CA 


VO LA 
LA O 
^ CM 


!>• tr- 
UA CM 
LA 


CA 

VO 


CA 


VD 

VO 

A 

A 

00 

LA 

VD 

CM 

LA 


LA 

VO 


LA 


LA 


o 

O 

LA 

LA 

LA 

A 

LA 

LA 




LA 

{>- 

CO 

in 

!>■ 

A 

tA 

VD 

CA 

LA ■ 

CO 



o 


CA 

{>■ 

O' 

^X) 

— 

O 

LA 

O 


vM 

A 

CO 

CO 

cv 

LA 

LA 

IC\ 



I> 


00 


CA 


CM 


K ci3 




contd 



65 



CM 

CO 


m 

£> 


VC 

u 

1 o 

o 

O 

!>■ 

in 

{>- 



T- 

- £>- 

CTi 

o 

CM 

l>- 

r- 

l> 

03 


KV 

K'x 

Ch 

cn 

CM 

CM 




ji 

■* 





CM 

CM 

1 ^ 

: "M 

■J CO 

o 

^r^ 


'M' 

in 

cn 

E>“ 

n CM 

m 

C'\ 

t<^\ 

in 

LTV 

CM 

?n 

0 



CD 

!>- 



CM 

CM 


N CM 

m 


Ch 

CO 

o 

CM 


K' 

C CM 

CM 

tn 

m 


in 

tn 

■vr 

-r- 


-r- 

£> 

f>- 



CM 

CM 











s {>. 

cn 

'M- 


m 

'M' 

CM 


CM 



CM 

^‘■ 

'■'M' 


CO 

{> 

' - 



cn 

cn 



CM 

CM 


CO 

o 

o 

CO 

Ch 


m 

m 

. ^ 


in 


Ln 

'M" 

Lf' 


{>- 




O^. 

C'v 



CM 

C>J 

' — 


















O 

in 

CM 

irc 

o 


CM 

-M- 

CM 

■r— 


CM 

CO 


r> 

[> 


cn 


cv 

CM 

cn 

cn 



rn 

rn 

f I 

^ i 

5 } 









f s 

CO 

03 

CM 

O". 

i> 

i>- 

VO 

inv. 

Ch 

cn 

O 

cn 

03 


cn 




C\i 

tr\ 

CO 

03 



Ln 

Ln 


CM 

CM 


cn 


CO 

cn 

O 

CO 

O 

f— 

r- 

CO 

VO 

00 

{>* 

o 


in 

in 

[> 

I>~ 



!>• 

GO 

i ^ 


•M” 

KV 

tn 


£> 

OJ 

o 

!>- 

CO 

O 

[>■ 

CM 

cn\ 

m 

o 

!>- 


^n 


CO 

i>- 

cn 

KV 


O 


1— 

•» 





CM 

CM 


o ■ 

o 

tn 

CM 

in 


Ch 


, — 

CO 

»— 

t-- 

CO 

O 

cn 

in 


CO 

o 


CO 

[> 

in 

'M- 

CTv 

O 


tn 

tn 





’M' 

in 


o 


CM 

!> 


03 

CM 


N 

in 

VO 

o 

CM 

Kj- 

in 

o 

.. — 

LTV 

tn 

cn 



vo 

VO 

«*— 

tn 


m 

KC 





VO 

VO 


VO 

GD CM 
CM K\ 
VO VO 


VO 

MO VO 


(T* O 

‘MD i> 


CO Cr\ 


CT\ KV 
‘vO 

CX) — 
(M 


03 O 


O -"'i” 
Lf\ Lf>. 


[>- 

m Lo 


o o 

'M- (T\ 

^ CM 
O O 
CM CM 


‘M- Ln 

v-j- 


0^ CM 
^ ir\ 


D- ^ 

'-02 


ch cn o cr\ cj^ 

J>- 'MD D- r- i>- 


'•M~ CO 


cn .- 

LO 'vO 


LO O 
CM TA 
Cvi CM 


CO 

^ C\j 
CD 


O CO 
‘"i- O 

'M- LPv 


CM IS^ 
VO VO 


'M' 

in 

cn 

o 

ir~ 

1— 

CM 

rn 

'r— 

T 

1 — 

t — 


in 

Cel 

03 

O 

,-yv 

U>/ 

03 

in 

m 



CM 

CM 




CO 

in 

in 


m 

VO 

in 

fr\ 




00 


OJ 

CO 

CO 


m-, 

cn 

CM 

CM 

m 

m 


-r- 

CM 

CM 


in 

in 

00 

in 

in 

ire 

n- 

in 

— 

CM 

o 

CM 

^n 

rn\ 

m 

in 


CO 

cn 

IC- 

m 

«. — 

, — 

m 

in 

in 

CO 

00 


*'4- 

■4" 

CO 

CO 





in 

VO 

cn 


CO 

VO 

CO 

m 

!>• 

in 

G\ 

i>- 

CM 

CM 

x— N, 

— s 

cn 

in- 

c 

n 

c- 

Ln 

o 

-r— 

O 

n- 

n- 

tn 

cn- 

VO 

• 

"M- 

in 

VO 

in 

in 

*I \ 

in 

'O 


.n 

CJ 

t'- 

'O 

m 

■r— 

• 


[ in 

in 





' 4 " 




CM 

•M 

t 

“M” 



! 














43 


1 














h 

m 

1 CO 

m- 


O 



o 


m 

-M" 

CO 

cn 

CO 

CO 

o 

i CO 

co 

o 

00 

CO 

n- 

VO 

GO 

!>- 

00 

CM 

\ 

m 

VO 

o 

1 J>- 


CO 

m 

•«— 

T— 

C\i 

CM 

in 

in 

in 


-.M" 




[ 







x~ 







CM 


i 














•r— 

! CM 















♦ ? 

1 K 

Ct 5 

rr| 

'cb 

w 

Ztj 

j-^i 

l~^ 

c 3 

rd 

o 

l- ! 

; -H 




""t j 

' O 
1 

1 













0 


1 














rH 
















fd 

03 

1 ^ 

vM* 


in 


VO 


£> 


CO 


cn 


o 



• 













CM 



contd, . 



1 







LA 

CA 

I>~ 


o 

LA 

LA 



LA 

CO 

00 

CM 

A- 



l>“ 

O 

!>- 

VD 


A“ 



KO 

T— 

T— 

CD 

O 


! 

\ 

LA 

\o 

vjD 


-r- 


1 ^ 

Ch 

LA 

{>“ 

LA 



1 ^ 




A~ 

CA 

i 

1 - 

1 


CM 

CM 



1 

KN 

CM 

LA 

LA 

A- 

MO 

tA\ 



CD 

O 

O 

lo- 

CA 

w- 

rs 

LA 

LA 

MO 




j 






CNJ 

-1— 

O 

00 

O 

co 

LA 

— , 


LA 

00 

-« — 

CA 

CA 

— ^ 

CM 

CM 

IA\ 





CM 

0'^ 

A- 

CO 

A- 

A“ 


CCi 

{A^ 


CM 

sj~ 

*0). 

,5 

0'\ 

•r— 

A- 

O 

tr^ 

LAi 

-M" 

CO 

» -Vf 

^C'^ 

LTx 

MO 

A~ 



03 

CM 


CM LA 

O 

LA 

LA •«“ 


o 


LA 

A~ 



CO 


AJ 

•nj* 





^ CM 

O 


£> CO 


o 


LA. 

LA 

CM 

lA 

lA 

CD 


LA CA 

o 

LC"‘. 

’M" A! 

LA 

O 


o 

LA 

A LA 

CD 

CO 


64 


V— : 

oo 

LA 

LA 

LA 



CO 

O 

lA. 

•«— 

o 




MO 

-r- 

LA 


00 

00 




<r— 


AJ 

A 


K-i- 

LA 

LA 

1 

A 

A 

LA 

CA\ 






▼— 

** 


OJ 



j 
















o 

O 

MO 

O 

CD 


•r— 


o 

A 



OO 

CM 


-t 


■ MO 

MO 

CO 

0- 


D- 

!A 

'S — 

** — 

o 

1 — 

MO 

'x'O 

MO 

0- 


1 LA 

1 

lA 

LA 

LA 









T— 




i 

00 

0 

A 


LA 



A 

A 

!>- 


o 

LA 


CA 

O 

A 

co 

O'^ 


lA 


MO 

00 

O 


CA 

o 

00 

o 


^ CD 

CD 

CA 

CA 


I ^ 

A 






00 









O 








X— 








•H 
















-P 












CA 

X— 

LA 


cd 

A 

LA 

A- 

o 

£> 


A 

MO 

CD 

CA 

— , 


LA 

00 

LA 

-H 

LA 

A 

A 

LA 



O 

00 

LA 

LA 

CD 


00 

CD 

CD 


X— 

xr- 





K's 

A 

X— 

X— 

\ 

•r— 

X — 

^ — 

*r— 

o 
















Q 




























A 

X— 

00 

LA 

U 

a^ 

co 

A 

X— 

LA 

A 

GO 

O 

LA 

A- 

N 

MO 

00 

CA 

A 

CO 

D- 

.p— 

LA 

(A 

X— 

LA 

CA 

LA 

LA 

00 


00 

A 

A 

O 


o 

X— 


LA 



A- 

O 

X— 



lA- 

00 

A 

O 


X— 

X— 





X— 

OJ 

X— 

▼— 




X— 

■c— 

Oj 
















-Pi 
















col 












LA 

!>- 

OA 

o 


A 

00 

X— 

X— 

LA 

CO 

CA 


A- 

CO 


00 

A 

l> 

A 


A 

A 


00 

LA 

o 

A 

LA 

X— 

CA 

CO 

CO 

D- 

00 

CA 


A 

A 

00 

X— 

A- 

LA\ 

A- 

LA 

CA 

A 


LA 

I>- 

CO 

D- 


"=M~ 

LA 


A 

X— 

A 

CO 

CA 

Kj- 

CO 





v— 













CO 

co 

A- 

,, 


00 

A- 

o 

o 

lA 

A 

A- 

00 

CA 

00 


LA 

LA 

o 

X— 


CA 

A 

A- 

A 

o 

LA 

CA 

CO 

O 

CA 

LA 

x 

A 

o 

CA 


CO 

X— 

LA 

A 

CA 

lA- 

LA 

(A 

X— 

A 


CA 

CA 

o- 

A- 


LA 


-r- 


A 

LA 


LA 


LA 




X 

X 














i 

lA X- 

X- A 

CO 

CD 

A- 


A 

LA 

CO 


CM 

U\ 


1^1 


A~ GO 

CA 

00 


LA 

MO 

A 


CO 

O 


. 

!M-| 

CO 

LA O..- 

X— 


CA 

O 


A 

nr> 

lA 

CA 

o 

* ' : 


LA CD 

[ 

CA CA 




X— 


X— 

A 

Cm! 


X— 

i 

\ ! 

1 












-P i 

i i 














y sJ 

1 A- A- 

CO -- 

UA 

^ — 

A- 

o 

A- 

CO 

CA 

MO 

X— 

' o 

o 1 


X- A 

■ A> o 

K^ 

MO 

A- 

A- 

LA 

LA 

MO 

MO 

Lr\ 


O 

8 

o 

o 

A 

A- 

A- 


<AJ 

X— 

LA 

CM 

A 

lA 

■K-j- 

^ 1 

1 


X— 





X— “ 






04 ! 


i 












x— 

<* i 

cm: 


Ixi' C0 

[xj 

rH . 

ixl 

CD 

!--J' 

!- H 

CJj 

(-T~l 
!•— J 

rj:s 

hj-i 

rH 

^ j 














0 : 














H ■ i 

1 



























a1 

P 

■ ■ 

A 

X— 


A 


LA 


LA 

EH 


A 

A 












contd 



Table 4.12 (contd. 


65 


LA 

O 

LA 

^A 



LA 

-r— 


NA 

LA 



CX) 


CD ^ 

o <y\ 

IT- CT\ 


VD KA 
cr\ VO 
C\J VO 
CM CM 


CM 

CM CM 

00 CO 


Lr\ -r- 

LTV. 

CM CM 


Lr\ 

GO 

VD ■ -<t 
<T\ 


AO VO 


UA 


^ {> 


^ {> 
CM 


O VO 


o 

CM CM 




Ln in !>- 00 


GD 


CM CM 


VO !>- 


'T ^ 


'X) CM 
Af- 


tSA 


CA -Ah 
CM COx 


(T\ VO 

00 00 


CM Ah 


^ 00 

CM 


LA 

VO O 


MD CO 
VO f> 


VO 


00 00 


Cv! -M" 
tA ?yA 


CD !>• 


-- l>- 

IPs ''vi' 


CT\ CM 
r- i> 


'MO 00 
■^4- 


fCAi 


00 




;MO 




l-<t\ 


;CSA 


eI 



CD 

— 


LA 

LA 

56 

!>- 

00 

OA 

CO 

M- 

co 

LA 

CM 


77 

77 

36 

'4- 

LA 

44 

Ah 

LA 

159 

LA 

LA 

VO 

CT", 

^A 

CM 

'vl" 

o 

LA 

577 

607 

CA 

O 

VO 

97 

CM 

VO 

LA 

120 

00 

'=4 

197 

O 

LA 

CM 

39 

CA 

CM 

O 

VO 

LA 

VO 

LA 

LA 

iA 

O 

00 

LA 

Ah 

VO 

VO 

ITS 

VO 

LA 

CM 

418 

1472 

671 

592 

956 

LA 

CA 

00 

00 

NA 

CM 

CO 

LA 

LA 

00 

CM 

1980 

4178 

4578 

LA 

CA 

LA 

00 

LA 

1623 

4958 

662 

o- 

o 

O- 

CM 

A» 

CA 

CM 

4685 

CA 

00 

LA 

CM 

CM 

O 

CA 

LA 

3541 

9599 

9617 

00 

CM 

O 

Lr\ 

00 

LA 


00 Ah 
Ah VD 
O VO 
fJA 


O CA 
VO -~ 
in 00 
CM 


A^ ‘sO 
^ CA 
CM CM 


LTV ^ 
VO VO 
LA LA 


LA VO 
K\ LA 
CM CM 


O GO 
O 00 


— 00 

^A VO 


CA <T\ 
LA T- 
03 CA 


w O 
O 00 

cn o 


VO 

O 

VO 

O 

CA 

VO 


CA 

CM 

-I 

LA 

LA 

LA 

Ah 

LA 



CA 

O 

CM 

'4" 

VD 


O 

.p— 


LA 

LA 

O 

O 

LA 

LA 


CM 

CD 

VO 

00 

^ 


CA 

CD 

LA 

00 

CA 

VO 

VD 

VD 

VO 

JA- 

LA 




LA 

LA 

00 

00 

Ah 


VD 

Ah 

lA 

lA 

VD 

£>“ 

lA 

CM 

o 


CM 

LA 

O 

-4- 

03 

OA 




CM 


CM 

!xi 

Cti 

!-7-i 

i'-’l 

■J3 

Fp} 

FH 



VO LA 
CA VO 
KN CM 
CM ^A 


i'-- £>- 

VO O 

LA. C- 


VD 


CD 


CA 


CM 


contd 



Table 4.12 ( con-td. 


66 




o 


in oj 
co\ 




c\J 


si 


oi 


OD 


M 




in 


^r^i 


CM 


CO 


^ O 

-‘vf 


MD 

ir\ 


CM CM 

tr\ tr\ 


^ LPv 
^ !S-\ 


CM 00 


CM 00 

Cn 

ITS ro 


tA CTi 
LPv MO 


o 

UP CM 
C7> ^ 


■M- CM 

IT- VO 
'vO UP 


VO .r- 

O CM 


W 


c<p 


CP LPv KP CM 

KP CM XT\ 

^ O 00 VO 

CM ^fP ^ ^ 


CM KP 


CP 
^^^ CM 


LCP UP 


!>- UP 


O 

KP ^P^ 


C 0 ‘ — 

LPv VO 


'M" OP 
O CP 


00 MO 
CO OP 


o- 

CP 


O 

KP O 
■ ^ 


VO 

O KP 
CP I'P 
V- CM 


KP rp- 

CP i>- 
^ LP 


CM p- 
CM {>- 


VO 

O UP 
"Vh KP 


KP CM 
CM O 
CM CM 


CM CO 
P- VO 


CO Pi 
^ !>. 
CM CM 


CM IPs 
UP DP\ 
CM CM 


VO 00 
00 CP 


O CM 
CP 

CM KP 


'vO KP 
^P MO 
UP 


CM CP 
CM 

KP KP 


CP ^ 

OO 00 


00 

I>“ VO 
K' I:P^ 


CP 
^p KP 
LIP 


VO CO 
00 CP 
KP 
LfP VO 


UP -M- 

bP ■»— 
UP 


O O 

T- KP 

UP VO 


LP\ 


o 

o 

VO 

IP- 

O CO 

,, 

00 



VO 

UP 

CM ^ 

CM 


K^a 

KP 



KP PJ 

1 >- 

!P- 








VO 

VO 

00 


IP- CP 

UP 

00 



UP 

I>- 




VO 

VO 

o 

IP- 

8 

6 

C^i 




K^ 

CM 



■» 

r- 

OO 

KP 

KP 

LPv KP 

KP 

CM 



LPn 

LPv 


CM 

04 

o 

o 

MD 

o 

CM 

PJ 

to 




CM 

KP rp 

LPv 

UP 

CP 

00 

CM 


"vP OJ 

OP 

VO 

i>~ 


UP 

VO 

'vp "'M" 


KP 






.r— 

CM 

IP- 

00 

CP 


CP CM 

liP 


MO 

VO 

o 

o 

O VD 


UP 



'M- 

VO 

KP 

CO 

CP 

o 

LPv 

CM 

KP 

O KP 

CO 

vq 

ip- 


UP 

CM 

O 

CM 


T— 


(P- 

I>- 

IPi LfP 

CM 

CM 






.r- 


KP 

CP 

VO 

l> 

00 UP 

EP- 


LPv 

t — 

VO 

O 

CP o 

vo 

CM 

CM 

CM 

tp\ 

IP- 

CP CP 

KP 




KP 

KP 


KP 

'M- 

OP 

r- 


O 

O -P- 

IP- 

VO 

MO 

VO 

00 

o 

CM CP 

00 

iP- 

CM 

CM 

D<p 


IP- UP 

IP- 

vo 



CM 

CM 




IP- 

VO 

CO 


CO UP 

UP 


IP- 

MO 


'r— 

VD !>- 

CP 

o 



CP 

P 

Pvt Cvj 

•* 

CM 



HH 

i- 1 -^ 

CO 

m: ei . 


cb 

VO 


IP- 


a> 

CP 






T— 

•v— 



contd 



67 


I 


Lr\ 

Ch 


Ch 

■r— 

VA 



CTi 

00 



00 

Ch 


Ln 

00 

■r— 



Ch 


£> 

UD 

LA 

OvJ 

Oh 

LA 

CV 








*r— 

O 

^~ 

CNJ 

C\J 

-r— 

CA 



DA 

KO 

UD 

o 

O 





OJ 

CM 

VO 

o 

LA 

LA 


<M 

CM 

o 

o 

O 

O 




LA 

CM 

VO 

CM 





-« — 


CM 

Cvj 







-i — 

CO 

•« — 

Ch 

VO 

LA 

LA 


CM 

CM 

o 

Ch 

!>■ 

O 







Ai 








-« — 







•r- 

E> 

LA 

<r- 

o 

Oh 

VO 


X— 

•«-“ 

•nO 

VO 

Ch 

Oh 








o. 

Ch 

VO 

VO 

Ch 

o 

VO 


CM 

CM 

o 

Ch 

I>- 

LA 




■r— 


T— 

„ 



Ch 

CM 

CM 

Ch 


Ch 

m 


Ch 


o 

VO 




f— 

CM 

LA 

LA 



CM 

Ch 

MO 


^ ! 

00 

o 

00 


LA 

Ch 

i>- i 


!Sh 

LA 

, — 

o 

CM 

Ch 1 





CM 

•r— 

1 




IA\ 

CM 

!>- 

MO 

o 

1 

00 i 


— N 

Ch 

{>“ 

MO 

CM 


VO 1 


[>- 

00 

£>- 

-<r— 

Ch 







«r— 

00 

00 

Ch i 








[ 



LA 

O 

LA 


00 

Ov] i 


N 

lA 

LA 

VO 

CM 

,, — 

c- { 


VO 

Oh 

Ch 


o 

'■M" 

LA { 




T— 


LA 

CO 

00 j 



Ch 

LA 

00 

CM 


; 



CM 

O 

LA 

'p— 

o 

00 


LA 

VO 

LA 

-« — 

CM 

VO 

O 



LA 

LA 

A- 

Ch 

Ch 

o 








,J — 



O^'i 

LA 

j- 

!>- 

MO 

o 

• 


LA 

CO 

o 

Ch 

CO 

A 

m 


CO 

VD 

LA 

^A 

Lxh 

00 

♦ 


C] 

CM 

LA 


LA 

A 









-¥=> 
















O 


lA 

!>- 

VO 

{>- 

00 

D- 

o 

LA 

CO 

T— 

MO 


i>- 




CM 

LA 

A- 

O 

<>— 

! 






-t — 

•< — 

Ai 

CM 








'r— 








• 









CM 

!i:l 

es 

|x| 

C±5 

Ixl 

cb 


' — -« 







(1) 








H 


j 






r-Q 








ci5 


! O 


•»— 


CM 




! CM 


CM 


CM 




68 


Table 4.13 Statistical properties of nistorical 
and generated flows for the entire 
Basin B3 

(Values in 10 M.Cft.) 



Maximum 

flow 

Minimum 

flow 

Mean 

Standard 

deviation 

del 
dable i 

Recorded and 
reconstituted 
flows for 

10 years 

576318 

217326 

374988 

126706 

276680 

Generated flows 
for 100 years 

993460 

119717 

383173 

166803 

269333| 


CiiAPTSa 
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RESULTS AIJU DISCU33ICU3 

I 

5 . 1 Introduction 

5.1.1 Gonoral 

Optimiza.tion and simulation techniques capahle of analy: 
in detail large water resources system are very valuable to ; 
planning engineer. The real complex system is represented b; 
model and the above techniques are applied to it. %il6 mod 
cannot make decisions, they can provide valuable informa-tion 
concerning the. construction and operation of a proposed set 
water resources projects. 

5*1.2 Itreamflow generation 

Historic records of streamflow are too short to includel 
all possible patterns of droughts and floods. The generatit 
of synthetic str-eamflows provides longer sequences which i 
permit more extensive analj^sis of the consequences of low anj 
high flows on the output of w.ater resources systems than do 
techniques which use only the historical record. In this 
study, a, monthly strea-mflox^ generation model developed by 
Hydrologic Engineering Centre is used to generate streamflowl 
The inf loxf data used varied from 7 years to 58 years for bas| 
BA and 7 to 1U x'-ea^rs for basin BB. On the basis of this infi 
mation, 138 years of flow has been generated for basin BA ani 
110 years for basin BB . The program takes about 50K memory 
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the DEC 1090 ayat-em. The CPU time for streamflow generation ; 
for basin BE. was about 2 minutes while for basin B3, which 
has a. larger system, it was about 5 minutes. The generated 
flows compare statistically wel.1 with the historical data 
as seen from Tables 4.10,4.13. The method adopted here, 
thus, can be used to get the data base for simulation studies, 

5.1.3 Optimization 

The genera.ted streamflows are used in the optimization 
study. The optimization technique used uses the out-of -kiltei 
algorithm which is one of the most efficient ways of solving ^ 
a network flow model. The program davoi.ope-d by the Texas 
Water Development Board is used with suitable modification. 

105 years of generated inflows are used. 7 runs were 
taken in the DEG 1090 computer system, using 15 years for 
each run. The GPU time for each run is about 14 minutes and 
the model roopuircs about 70K memory. The results of the 
simulation presented in this report include the demand 
deficits at each node, the surplus water at spill nodes, annual 
powex* genera.ted ana end of the month stojeages. The frequency ; 
analysis of the demand deficits give a method of evaluation : 
of individual projects with respect to their expected perfor- 
mance. The statistical analysis of the spill.s indicate the 
quantity of water available after satisfying the in-basin 
demands. The end of the month storages are used to develop 
rigid operating rules for the reservoirs. 
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5*2 Results of Studies and Conclusions 


In Indi 

.a, 

irrigation projects ar 

c planned such 

that 

demands are 

to 

be met with 75 percent 

d ep end ab i li ty 

1 P 

100 percent 

of 

the domand should be in 

et at least 75 

percent 

of the time. 


The frequency anaj.ysis 

o:' dSTnaiidg met 

or 


exceeded is shown in Table 5.1. Table va7Lues indicate 
the percentage of time the stated frequencies of demand 
being met or exceeded. It is found from the tab]e that 
nodes 4 and S have large deficits. I'or node 4, only 
30-48 percent of time, 95 percent of demand is met or 
exceeaod and for node 6, only 39-05 percent of time, 95 per- 
cent of demand is met or exceeded- This shows that the 
two projects represented by those two nodes have not been 
properly planned- The demand to be met is much higher than 
the inflow a.vailable at tnesc two sites. Por all other 
resorvoirs the design cr’iterion has been met or nearly met. 
lig. 5-1 shows a plot of parcentage of demand met or exceeded 
with tiie percentage of time for two selected nodes. ^ 

Table 5-2 shows the statistical, properties of spills. 

The amount of .a.ter availa.ble for out of ba.sin transfer at 
nodes 1 and 17 after meeting the in basin demands on a monthly 
basis is indicated in the table. At node 1, it is found that 
surplus flows occur only during monsoon months while at 
node 17, a small quantity ox water flows out even during nony 
monsoon months. This is due to, water coming from the sub- I 

basins in between nodes 1 and 17- As seen from the table, I 
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Fig. 5*1 Frequency analysis of demand met at selected nodes 
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Note ; 'rhe demands at nodes 10 and 12 are zero. j?iode 20 is an artificial spill node. 

faPle values indicate the percentage of time the stated frequencies of demand 
oeing met or exceeded. 



Table 5.2 Statistical properties of spills 
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an aver-age annual surplus of 273.82 TMG is available at 
node 1 for out of basin transfer. 

ilnother output of the optimization model is the end of 
the month storages for all the reservoirs in the system. 

These are classified for the three hydrological stat^ namely 
dry, average and net. The mean value calculated for each 
month and each hydrological state form the rigid rules for 
reservoir operation. Table 5.3 shows the.se rigid rules for 
all the reservoirs in the system. It is seen from the table 
that the dry year operating rule during monsoon period is 
higher than that of average and wet years. This indicates 
the carryover storage brought from the end of the previous 
year, during non-monsoon period, the dry year rule is at 
t;ie lowest indicating that the reservoir storage has to be 
lowered to meet the demands. The wet j'-car rule during mon- 
soon is the lowest indicating very small quantity of carry- 
over storage, during non-monsoon, it is the highest indicating 
storage preserved for future use. These trends are observed 
in most of the reservoir opera.ting rul.u developed here. Plots 
of rigid operating rul.e for node 1 and 4 are shovm in 
Pig . 5 . 2 ( a. ) and ( b ) . 

The results of the study also contain the annual pov/er 
geiioratod for each of the povfer stations in the system. The 
annual poxirer generated is statistically analysed to find the 
mean and standard deviation and the resul-ts are presented in 


Table 5-4. 
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Fig. 5-2 (a) Rigid operating rule for Node 1 








Table 3.4 Statistics of annual power generated 
(all values are in 10 kilowatts) 


S' ode 

Power 
Ifo . 

Sta.tion 
link No. 

Mean 

Std. lev. 

1 


22 

24447.30 

2486.62 

12 


23 

1689.71 

588.69 

21 


24 

765.62 

48.55 

22 


25 

2864.93 

540.02 

23 


26 

30962.54 

1400.02 

24 


27 

654.51 

28.40 

25 


28 

15896.06 

998.26 

26 


29 

947.21 

87.63 

10 


30 

870.56 

523.49 

1 1 


31 

372.23 

193.74 

7 


32 

4551.73 

2183.90 



Ill coiiclusiorij it may be stated that large water 
rosoixrces system can be analyzed by using the technique 
followed in this study. The reliability of the results of 
such studies depend not only ^on the proper representation of 
the physical system by a model;, but also on the use of 
reliable data. The model requires to be updated periodically 
to reflect the chajiges that take place in the basin. 
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